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Abstract 

We calculate the heat kernel coefficients for the electromagnetic 
field in the background of a dielectric cylinder with non equal speeds 
of light inside and outside. The coefficient 02 whose vanishing makes 
the vacuum energy of a massless field unique, turns out to be zero 
in dilute order, i.e., in order (e — 1)^, and nonzero beyond. As a 
consequence, the vanishing of the vacuum energy in the presence of 
a dielectric cylinder found by Casimir-Polder summation must take 
place irrespectively of the methods by which it might be calculated. 



1 Introduction 

The vacuum energy of the electromagnetic field in the presence of a dielectric 
body is a topic of considerable interest since Schwingers attempt to explain 
sonoluminescence. Although this seems unlikely to be possible, mainly due 
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to different time scales involved, open questions connected with the renor- 
malization remain. 

The divergent part of the vacuum energy can be expressed in terms of 
the corresponding heat kernel coefficients. By means of the renormalization 
procedure the divergent part must be subtracted from the ground state en- 
ergy and added as counter term to the classical energy of the background. 
This general procedure does not work straightforwardly for a dielectric body. 
The point is that there is no classical energy which may be associated with 
it. Note that there are no electric or magnetic fields present, instead there 
is merely the ability of the body to get polarized upon application of an 
electromagnetic field. One may overcome this problem by introducing for- 
mally the necessary terms for a classical energy and put zero the coefficients 
in front of them. This is similar to semiclassical gravity where one assumes 
the coefficients in front of the terms quadratic in curvature to be sufficiently 
small. But there is one more problem associated with the vacuum energy of 
a dielectric body. As the field is massless, there is no normalization condition 
to get a unique ground state energy as long as the heat kernel coefficient 02 
is nonzero because it comes with the logarithmic divergence, see ^ for a 
detailed explanation. 

In [U the heat kernel coefficients for a generic smooth dielectric back- 
ground and for a dielectric ball have been calculated. The coefficient 02 
turned out to be nonzero in general. For small deviation of the dielectricity e 
from unity (resp. from its value at spatial infinity) 02 is zero to first order in 
{e — 1) and nonzero starting from second order, [e — l)^, which is called dilute 
approximation. For the dielectric ball it is zero in second order in addition. 
In this way, for a more singular background configuration {e{r) = 6Q{R — r) 
has a jump) the singularity is slightly weaker. It is this zero of 02 to order 
[e — 1)^ which allowed to get a unique result for the vacuum energy of the 
dielectric ball in dilute approximation, 

which was obtained in by Casimir- Polder summation in 0, by a perturbative 
setup in and by mode summation in 

In the present paper we consider a dielectric cylinder and calculate the 
heat kernel coefficients a„ for n < 2. This problem is technically more 
involved than the corresponding calculation for the ball. Although the Jost 
function of the corresponding scattering problem which we use is expressed 
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in terms of Bessel functions the photon polarizations do not separate and 
there is an additional dependence on the axial momentum. 

The vacuum energy of a dielectric cylinder had been calculated in 
in dilute approximation by summing the Casimir-Polder potentials and it 
turned out to be zero, 

£f = 0{is- 1)3) . (2) 

In the present paper we show that the corresponding 02 is zero in this (dilute) 
approximation and conclude that the vacuum energy in dilute order is zero 
irrespectively of the method by which it is calculated. 
In this paper we use units with h = c = 1. 

2 Basic formulas 

In order to regularize the vacuum energy of the electromagnetic field we use 
the zeta functional regularization and express the energy in terms of the 
corresponding zeta function, 

W = ^^c(.-i). (3) 

In the presence of the dielectric cylinder the spectrum of the electromagnetic 
field is completely continuous and we use the representation of the regularized 
ground state energy derived in for a spherically symmetric background 
which was in |^ rewritten for a cylindrically symmetric background. After 
dropping the Minkowski space contribution we obtain 

c(.) = ^f: J^JdkiP-kir^inMk,h). (4) 
'=-~ 

Here, Ai{k, k^) is the Jost function for the scattering of electromagnetic waves 
off the cylinder with radius R on the imaginary axis. It is given by 

MKh) = -^{Aj'\k,k,)Ar{k,k,) (5) 
+/2 - kl) kl (1 - ^2)2^-2 [i^^qR)Ki{kR)f] 

with 

A™ (A;, k,) = fiikl'iiqR) Ki{kR) - ii^qlMR) K[{kR), 
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Aj^{k,h) = e^kI[{qR)Ki{kR)-e2qh{qR)K[{kR), 
Ar{k,k,) = ie2(^-^)gfc(£iA; + £2g)(/iiA; + /i2g), (6) 



and the notations q = ya'^k'^ + k^l — a^), a = C2/C1; and the speeds of 
hght inside and outside are q = l/^/eiJIi {i = 1,2). This formula can be 
found in textbooks on classical electrodynamics, see also @]. 

The heat kernel coefficients are known to be defined by the residues of 
the zeta function at the corresponding points s, 

a„ = (47r)3/2Res,=|_,r(.)C(s). (7) 

These residua result from large all, the momentum k, the axial momentum k^ 
and the orbital momentum /. To factorize the orbital momentum dependence 
for the terms with / 7^ we use the uniform asymptotic expansion of the 
Bessel functions Q . Being inserted into Eq. it gives raise to the expansion 

In Mlk,lk,)r. (8) 

j=-l,0,l,... ' 

for fc, / — > 00. Substituting the latter into (^) previously changing the vari- 
ables k ^ Ik and kz — >■ Ik^ so that k and k^ were fixed, we split the zeta 
function into parts according to 

r(s) Ccyi{s) = A{s) + 2 + Ao{s) + Ai{s) + A2{s) + A^is) + ...}. 

(9) 

Here A{s) results from / = 



and below we will use the asymptotic expansion of the Bessel functions for 
large argument in it. The Ai{s)^s are generated by the expansion (§) 



00 , , 00 



= jr(i-.) J iri*<* - (11) 

For the heat kernel coefficients up to 02 it is sufficient to include i up to 3. 
The summation over I (from 1 to 00) was carried out in Ai{s) and resulted 
in the Riemann zeta function. 
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The expressions for the Di are simple for i = —1 and i = 0, but sufficiently 
involved for i = 1, 2, 3 to be banned into the appendix. Here we note 



D. 



, q l + ^/lTq^' 
m — 




4 P gSyi^^/l + p 



{eik + e2q){^.k + /i2g) 



(12) 



In the next section the functions Ai [i = —1, ..3) and A{s) are considered 
one after another in order to find their residues contributing to the heat- 
kernel coefficients. 



3 Calculation of the heat-kernel coefficients 



In this section we give a detailed calculation of the heat-kernel coefficient 02. 
The results of the analogous calculations for the junior coefficients a„ with 
n < 2 are listed in the end of the section. 

We start from A_i{s), where it is easy to carry out the integration ex- 
plicitly 



A_l(.) 



[a 



(13) 



1 - 2s 

We are interested in the residue of A_i{s) at the point s = —1/2. As A_i 
has no pole at this point, the residue equals zero. So, this term of the 
expansion (^ doesn't contribute to 02- 

The term Aq{s) is more complicated. To simplify the integration we make 
the change of variables k k, kz —>■ krj. It gives 

C(2s - 1) 



Ao{s) 



71 



cl'Til - s) 



± 

I' 



The function 
Doik,v) = 



In 



'1 + 7^ 
1 + F 



00 

Jdkk-^'+^ 


d 1] d 





dk k dr] 


.1/4 


( i + e v 


1 +/^27' 


^l+72pyl 



Doik,r]). 



1/41 
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+ 



In- 



with 7 

ity, 



[1 — a'^)r]'^, has the following asymptotics at zero and infin- 



Do{k,r])\k-^o = In 



(/xi + 7V2)(£i + 7'^2) + (1 - V') (1 - 



-Do(A;,?7)|jfc- 



7(£i + £27)(a*i + A*27) 

(l-Q^)V(l-77^) 

2c27^(/ii + Ai27)(^i + ^27) k 



Therefore the integral over k in Aq converges in the range — 1/2 < s < 1/2. 
Adding and subtracting from the integrand its asymptotics at /c — > 00 we 
obtain the analytic continuation into the vicinity of the point s — —1/2 



Aois) = 



C(2s - 1) 



TT crr(i-s) 



X 



r dk 



d r] d 
dk k dri 

d 



772(1 -r;2)(l-a2)2 



(14) 



-2. + l)+,- 



r(-2s + i)r(2s + 1 



2ci72(^i + /i27)(£i + £27) (A; + l)2 
r/2(l-r/2)(l-a2)2 | 



C27^(ah + At27)(^i + ^27) 



The first integral is now analytical around s = —1/2, the second one 
contains the pole. But this pole is canceled by the Riemann zeta function 
C(2s — 1) which equals zero at s = — 1/2. As a result we have 



res Ar^(s) = 0. 
Analyzing now the asymptotics of Di{s) at zero and at infinity, 

l.O ,N.O „, ©1 



(15) 



Di\k^0^^i^"-^)k' + O{k'), D, 



k 



+ o{k-'), 



where 



4 1 7 £1 + £27 A*i + A*27 
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one comes to the conclusion that the integral over k in 



Ms) 



C(2s) K 



71 cfn 



2 1 



oo 


" 9 


1] d 







k dr] 



Di{k,T]), 
(16) 

converges in the strip < s < 3/2. The analytic continuation to the point 
s = —1/2 is the following 



AAs) 



C(2s 



7rr(l - s) 

1 , oo 



dr] 



[ok kor] 



Di{k,r,) 



-T(s) 



drj 



2\s 



(17) 



Here the asymptotics of Di{k, 77) for k ^ 00 was added and subtracted. Both 
terms in (|l^) are analytic around s = —1/2. Therefore res^ Af"^(s) = 0. 

Similar arguments for ^2(5), where the integral with respect to k exists 
at -1/2 < s < 3/2, give 



res ^2(5) 



dr] 



' 32 7r3/2 J 



(18) 



X 



72 (/ii + /i2)2 

The next term in (|^) we have to consider is ^3(5). The integrals in ^43(3) 
converge at s = —1/2, while the Riemann zeta function has a pole at this 
point. Thus 



res A^is) 



TC 



3/2 



J dr] y^l-r/2 J dk k^ 



d r] d 
dk k dr] 



(19) 



Now we are left with the term A{s) in Eq. ([T0|) 
1 



A{s) 



K 



2s-l 



cfTTT l - S) 



,2\-s 



00 

Jdkk-^'+' 


" d 


r] d 





dk 


k dr] 



\nAo{k,r]). 

(20) 
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The asymptotic behavior of In Ao(A;, r^) at = and A; = oo is the following 
lnA(A:,r^)U=o = In , ^ ^ + 2 (1 - 7) A: + ^(A;^), 

1 \(U W 1 -gi7-3£i + 3£27^ + g27 1 

InA A;,!] fc=oo = . t (21) 

8 [ei + £2 7) 7 «^ 

^ -3 gj 7^ + 4 gi £2 7' + gj 7' + e\ 7^ - 1 



1672(51 + £27)' fc' 
^/ii} + 0(A;-=^). 

The integrals over k for the first an the second terms in the square brackets 
of (pOD exist correspondingly in the ranges < s < 1 and < s < 1/2. We 
divide the integration with respect to k in two parts /q°° = /g^ + /i°°, where 
the first one is regular at s = —1/2, and the second one diverges in the upper 
bound. Keeping in mind that only the singular part gives contribution to 
the residue we omit the term convergent at s = — 1/2 . To obtain the 
analytic continuation to the point s = —1/2 one should add and subtract 
from In A(A;, rj) two terms of its expansion in infinity. After that it is easy to 
find the residue 

' 

" 3/xi7^ - /XI72 - 47Vi/i2 + 3/^2 - ^2^2 ^ ^ 
2? ] ^2 r ^ gj 

which up to the multiplier —1/2 coincides with (|THp 

Summing up we see that only the functions A2, and A have nonzero 
residues in s = — |. But as the residues of A2 and A are mutually canceled in 
d^), the heat-kernel coefficient 02 is totally defined by the contribution from 
^3. The integration over rj and k in terms of elementary functions seems 
impossible in (|T^). For /Xj = 1 simplifications occur, but still we are able 
to analyze the behavior of the coefficient 02 as a function of the velocities if 
light only numerically (see Fig.l). 

In the limit of small differences of the velocities of light (dilute dielectric 
cylinder) one might expand res ^3(5) in powers of Ci — C2. The expansion 

starts from the third order 

res Ccyiis) = --^^(ci - C2f + .... (23) 



'-1/2^'y''"' 1408 0Fci 
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Figure 1: The heat kernel coefficient a2 for the dielectric cylinder as function 
of the ration Ci/c2 of the speeds of light inside vs. outside the cylinder 

The heat kernel coefficient 02 is 



0-2 



TT 9 



(Cl - C2f + . . . . 



(24) 



The heat-kernel coefficients with n<2 are obtained in the same way. For 
arbitrary velocities of light inside and outside the cylinder one gets 



03/2 



C2 



ai = 



2 7r=^/2i? {l-a'^y 



<A (a«i + A«2)(£i + £2)' 



27r 8 Jdk j rj d 

00^'' 



A (A;, 7;), 



+ 



drji] 

C2 J y/l-T]'^ 



or] 



47'^£2/i2 



(/Xl + /i2 7') (£1 + £2 7') + (1 - - C2-' 



3 TT^/^ e^/Xi - 4£i/ii£2Ai2 + £2/^2 + ^llA + ^2Ail 



8i? 



(£1 +£2)^ (a*! + A*2)^ 
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4 Conclusions 



In the preceding sections we calculated the heat kernel coefficients a„ with 
n < 2 of the electromagnetic field in the background of a dielectric cylinder 
for different speeds of light inside and outside. The coefficient a2 vanishes in 
the dilute order. 

The result of a non vanishing a2 beyond the dilute order, obtained in 
m for a dielectric ball and here by Eq. (^) for a cylinder implies an yet 
unsolved problem. On the one hand side the vacuum energy of the electro- 
magnetic field in the background of a dielectric body is a physical quantity 
and at least in principle measurable. On the other hand side even after re- 
moving the ultraviolet divergences by hand, it is not uniquely defined due to 
the arbitrariness resulting from the logarithmic divergence. The only expla- 
nation which seems meaningful states that the problem itself is ill defined. 
The setup of a material body characterized by a position dependent a 
step-like e{r) = eQ{R — r) for instance, is clearly an idealization. Another ex- 
ample for an idealization are conductor boundary conditions (plane mirrors 
or a infinitely thin spherical shell for example) delivering a uniquely defined 
vacuum energy. 

So this idealization does not work for a dielectric body. This is surprising 
in view of the whole area of electromagnetic phenomena like scattering of 
waves not connected with vacuum energy where it works well. In order to 
attack this problem one has probably to take some microscopic model for the 
dielectric body, a lattice of harmonic oscillators for instance and to calculate 
the vacuum energy in that background beyond the dilute approximation. 
Then one had to investigate a limiting procedure to turn to the continuous 
dielectric body. One should expect that a dependence of the vacuum energy 
on some parameter like lattice spacing or plasma frequency remains resulting 
in the corresponding logarithmic term connected with a2- There are recent 



attempts to do so, see [10, 11, 12 
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Appendix A 



In this appendix we give a list of functions entering the expansion (P). They 
are written in terms of the variables k and t] with making use of the notations 
ti = t2 = l/VTT¥, 7= Ja^ + {l-a^)r]^ and 



Diik,ri) 
D2{k,v) 



^1, 



(25) 



where 
^1 = 



1 {{\-tX){\-t\ 



4.2 f2 



Sljll 



tl) 



+(eiU2 + e2Ui) (-— + — - — + ^ I 
^2(1 



12 



^2 7 ^2 7 ^1 "I" 



+£2 /i2 



■t2 



12 



tl 5 
TT2 



+£2/i2(l-r?')^?'(l-a')'tit2 T7T^2-T + T -TTT^I 



12 



4 4 



12 



12 



(26) 



To 




tl (1 -ti) / 3 

'^^^t2(l-t?) U'^'^ 



1 tltl [-^^u,^l-tl) 

7 g 35 3 3 ^ 2 205 ^ 

— tl to -l- ti to ~h — to ~h tr 

48 ^ 144 ^ 2 32 ^ 288 ' 



+ 



48 

-(£l/i2 + £2/^1) 

35 
~ 288 

3 2 _ 203 6 _ 41 4 205 
^32 ^ ~ 288 ^ ~ 48 ^ ^ 288 

+e2^l2{l-rj^)rf{l-a'ftlt2 



3 2 
— ^2 

32 ^ 



^2 



iJ vl6 

^ 2 32 

1 i 3 

— t^ to -|- — t^ to ~h 

16 48 ^ 

48 ^ 288 ^ 32 

3 7 

tits tlto + 

48 ^ 

205 5 

^ 32 
1 



5 3 _ 41 4 

tl to to 

16 ^ 48 2 



vl6' 



13 
' 16 


''1 


203 

' 288 


— d + 

48 ^ 


1 s 
— t? 
48 ^ 


'0 






13 
— t 
16 


4 
2 


5 o 

— t? 
16 ^ 









'l6 



tlt2 



5 2 



32 
144 



35 3 3 

144 1 2 



11 4 

48 ^ 



11 



5 o 

5 2 

+ 32*' 



25 



3 3 O ' 



205 6 
288 ^ 



41 



205 



48 ^ 288 



''1 



(27) 



205 



/2 /2 



tl (1 - tj ) 



1015 6 3 13 4 65 4 . 

5 2 s 21 , 3671 
128 



zuo g 203 6 1015 
384^^^2+1152^^^'" 3456 



-tl) 



5 2 41 4 

128 2 64 2 



1435 o 

H t? to H 

3456 ^ 2 128 

4543 7 22715 g 35 

tl H to + 

1152 2 10368 2 ^ 

15 ^3 _ 469 ^7 1103 ^5 
"128 ^ ~ 128 ^ +"640" ^ 

+ (£l//2 + £2A^l) 

35 



,2 ,3 _ ^ 

''1 ''2 IOC 



— 4 

128 2 



.^5 

1920 2 
3 2 287 3 4 21385 g 

384 ^2 " 576 ^2 + 10368 



35 6 
1152 ^^^2 



4 



/ 3 2 11 4 '^<-' 6 

IT28^^^'"T92^^^2 + 

35 g 2 H 4 H 4 3 
tl to ~\~ tl to — — z — ti to 

3456 ^2 192 ^ 576 ^ 2 
1 .9 .s . 9 .s 293 



35 



1152^'^' 

35 3 ^2 

576 '2 ~ 3456 ^ ^ ~ 128 ^ ' 
2 3 y 3 293 5 789 7 3115 1 

128 ^2 + 128 ^2 ~ 640 + 1152 " 10368 ^ 

11 3 4 3115 9 93 " 

576 ^2 + 10368 ~ 128 ~ 



787 



4 ^^^4tl 



,0, 7 293 
128 ^ 1152 ^ . 

t2 (1 - 1?) 3 2 13 4 203 6 205 ^ 



t2(l-tD / 3 2^ 13 
'ii(l-ii) V 128 ' g4 
1435 6 3 41 4 287 4 3 1015 36 15 2 

3456 ^2 64 ^ 576 ^ 2 3456 ^ ' 128 ^ 
35 2 3 15 3 1103 5 469 7 21385 9 53. 

384 ^ 2+128 2 ^ "640" ' + 128 ' ~ 10368 ' + 128 ^ ^ 
65 ,^ ,4 22715 ,q , 21 ,3 , 4543 ,7 3671 ^5 ' 



64 

2 3 ^ 3 

^^^2+128^2 
_t3^4 _^^9 , 21 ^3 , 4543 
192 ^ 2 10368 ^ 



l7 



— t? H r-i 

128 ^ 1152 ^ 

+ (1_^2^^2(1_^2^2 



.3 +2 

'2 



■ a'^ye2lJ'2tit2 

25 6 3 41 . 
3456 ^ 2 



41 
192 

+3 



128^^^^ 



1920 
41 



ttt2 



205 6 1025 

t? t2 + 

1152 ^ 3456 i„_ 
5 2 25 2 3 21 3 3671 
128 ^ 384 ^2 128 2 



205 4 3 
576 ^^^2 



41 4 205 

t^ to ~\~ ^1 t- 

192 2 1152 



^5 

1920 2 



1025 6 

3456 ^ 2 
4543 7 22715 

1152 ^2 + 10368 



12 



25 3 2 205 3 4 22715 g 
3671 



1920 



21 
128 



4543 
1152 



(28) 



and the denominator is 
{l-tl){l-tl) 



f2 f2 



t2 (1 - ti) 



ti{l-£) , , t2 

^1/^1 rr^^ — 12T + (£iA«2 + £2A«i) +£2/^2 



ti 1 



t2 



(29) 
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